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O ' Abstract 

A consistent folding model analysis of the (AS = 0, AT = 1) charge exchange (p, n) reaction 

^ ; 

measured with 48 Ca, 90 Zr, 120 Sn and 208 Pb targets at the proton energies of 35 and 45 MeV is 
done within a two-channel coupling formalism. The nuclear ground state densities given by the 
Hartree-Fock-Bogoljubov formalism and the density dependent CDM3Y6 interaction were used 
as inputs for the folding calculation of the nucleon optical potential and (p, n) form factor. To 
have an accurate isospin dependence of the interaction, a complex isovector density dependence of 
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the CDM3Y6 interaction has been carefully calibrated against the microscopic Brueckner-Hatree- 
Fock calculation by Jeukenne, Lejeune and Mahaux before being used as folding input. Since the 
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\ isovector coupling was used to explicitly link the isovector part of the nucleon optical potential to 



the cross section of (p, n) reaction exciting the + isobaric analog states in 48 Sc, 90 Nb, 120 Sb and 
208 Bi, the newly parameterized isovector density dependence could be well tested in the folding 
model analysis of the (p, n) reaction. The isospin- and density dependent CDM3Y6 interaction 
was further used in the Hartree-Fock calculation of asymmetric nuclear matter, and a realistic 
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estimation of the nuclear symmetry energy has been made. 
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I. INTRODUCTION 



The charge exchange (p, n) reaction is well known as an effective tool to excite the isobaric 
analog of the target ground state. Such an isobaric analog state (IAS) has essentially the 
same structure as the ground state of the target except for the replacement of a neutron by 
a proton and, hence, differs in energy approximately by the Coulomb energy of the added 
proton. In the isospin representation, the two states are just two members of the isospin 
multiplet which differ only in the orientation of the isospin T. The similarity of the initial 
and final states of the (p, n) reaction makes this reaction very much like an elastic scattering 
in which the isospin of the incident proton is "flipped". Indeed, the isospin dependent 
part of the proton-nucleus optical potential (OP) was used by Satchler et al. jl| some 40 
years ago as the charge exchange form factor in their study of the (p, n) reaction within the 
distorted wave Born approximation (DWBA). In general, the central nucleon- nucleus OP 
can be written in terms of the isovector coupling [2J as 

U(R) = U (R)+4U 1 (R)^, (1) 

where t is the isospin of the incident nucleon and T is that of the target A. The second term 
of Eq. ([T]), known as the Lane potential, contributes to the elastic (p, p) and (n, n) scattering 
as well as to the charge exchange (p,n) reaction While the relative contribution by 
the Lane potential U± to the elastic nucleon scattering cross section is small and amounts 
only a few percent for a neutron-rich target , it determines entirely the (Fermi-type) 
AJ n = + transition strength of the (p, n) reaction to isobaric analog states. Therefore, the 
(p, n) reaction can be used as a reliable probe of the isospin dependence of the proton-nucleus 
OP. 

The nucleon-nucleus OP has been studied over the years and there are several "global" 
sets of the OP parameters deduced from the extensive optical model analyses of nucleon 
elastic scattering, like that by Becchetti and Greenlees Q], the CH89 global OP Q], and a 
recent systematics by Koning and Delaroche which covers a wide range of energies (from 
1 keV to 200 MeV) and target masses (24 < A < 209). Although parameterized in the 
empirical Woods-Saxon form, these global systematics are very valuable in predicting the 
nucleon-nucleus OP when elastic scattering data are not available or cannot be measured 
which is the case for the unstable, dripline nuclei. Given a large neutron excess in the 
unstable neutron-rich nuclei, it is important to know as accurate as possible the isospin 



dependence of the nucleon-nucleus OP before using it in various studies of nuclear reactions 
and nuclear astrophysics. We recall here the two main methods used so far to determine the 
isospin dependence of the nucleon-nucleus OP: 

(i) study the elastic scattering of proton and neutron from the same target and measured 
at about the same energy, where the isovector term of Eq. ([T]) has the same strength, but 
opposite signs for (p, p) and (n, n) elastic scattering; 

(ii) study the (isospin-flip) charge exchange (p, n) transition between isobaric analog 
states. 

Since there are not sufficient high-quality (p, n) data available for a wide range of target 
masses and proton energies, the empirical isospin dependence of the nucleon-nucleus OP 



has been deduced 
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based mainly on method (i). While these three global nucleon 
optical potentials have been widely used in predicting the nucleon-nucleus OP in numerous 
direct reaction analyses within the DWBA or coupled-channel (CC) formalism, their isospin 
dependence has been rarely used to study the charge exchange (p, n) transition between the 
IAS's. The (phenomenological) Lane potential U\ has been studied in details so far at some 
particular energies only, like the systematics for U\ deduced from IAS data of the (p, n) 
reaction measured at 22.8 9] and 35 MeV [lo|. Therefore, it is necessary to have a reliable 
microscopic prediction for U\ by the folding model, to reduce the uncertainty associated 
with the isospin dependence of the nucleon-nucleus OP. 

Another very interesting microscopic aspect of the Lane potential is that it provides a 
direct link between the isospin dependence of the in-medium nucleon-nucleon (NN) inter- 
action and the charge exchange (p, n) reaction, so that accurately measured (p, n) cross 



section can serve as a good probe of the isospin dependence of the NN interaction [TJJ if 
the wave functions of the involved nuclear states are known. Furthermore, within the frame 
of many-body calculation of nuclear matter (NM), the asymmetry of the NM equation of 
state (EOS) depends entirely on the density- and isospin dependence of the NN interaction 



12, ll3J]. This asymmetry is actually determined by the NM symmetry energy S(p) which is 



defined in terms of the NM binding energy B(p, 6) as 

B(p, 5) = B(p, 0) + S(p)5 2 + 0(5 A ) + ... (2) 

where 5 = (p n — p p ) / p is the neutron-proton asymmetry parameter. The contribution of 
0(5 A ) and higher-order terms in Eq. (T5]), i.e., the deviation from the parabolic law was proven 
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to be negligible [12 , Q. The knowledge about the nuclear EOS ([2} is well known to be vital 
or the understanding of the dynamics of supernova explosion and neutron star formation 
161 ]. The NM symmetry energy determined at the saturation density, E sjra = S(p ) 



1U1 tuc 

3,H, 



with p ~ 0.17 fm 3 , is widely known in the literature as the symmetry energy or symmetry 
coefficient. Although numerous nuclear many-body calculations have predicted E sym to be 



around 30 MeV \\2 



13|,H,ll8|, a direct experimental determination of E sym still remains 
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22 1 or fragmentation data 



a challenging task. One needs, therefore, to relate E, 
quantity like the neutron skin in neutron-rich nuclei 
of the heavy-ion (HI) collisions involving N ^ Z nuclei 23j, |2J, |25j. In our recent study 
of the IAS excitation in the jo( 6 He, 6 Li*)n reaction using the folded Lane potential U\ for 
the charge exchange form factor 26], we have shown how the NM symmetry energy can 
be linked to the charge exchange (p, n) transition strength and, hence, be probed in the 
folding model analysis of the (p, n) reaction. To extend the folding model study of the (p, n) 
reaction to heavier targets to validate the conclusion made in Ref. 26] for the NM symmetry 
energy, we have studied in the present work the quasi-elastic (p, n) scattering measured by 
the MSU group for 48 Ca, 90 Zr, 120 Sn, and 208 Pb targets at the incident proton energies of 
35 and 45 MeV [ll|]. For a detailed probe of the isospin dependence of the in-medium A^A" 
interaction, a (complex) isospin- and density dependence of the CDM3Y6 interaction 27] 
has been carefully parameterized based on the Brueckner-Hatree-Fock (BHF) calculation 



of nuclear matter by Jeukenne, Lejeune and Mahaux 28j. While the isovector part of 



the nucleon-nucleus OP in the NM limit has been investigated in numerous BHF studies 
(see, e.g., Ref. |29J and references therein), the isospin dependence predicted by such BHF 
calculations was rarely tested in the DWBA or CC analysis of the charge exchange reaction 
to isobaric analog states. Our present folding model study provides, therefore, an important 
method to link the BHF results to the descriptions of the quasi-elastic (p, n) reaction. By 
using the Lane potential U\ to construct the charge exchange (p, n) form factor based on 
the isospin coupling, it is also straightforward to probe the isospin dependence of existing 
global nucleon-nucleus OP 

flag 

. In the present work, the description of the considered 



(p, n) reactions by the three global nucleon optical potentials 



flfl 



has been given, with 



a detailed comparison between the results given by the CH89 potential [7] and those of the 
folding model analysis. 
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II. IAS EXCITATION, LANE POTENTIAL AND ISOSPIN COUPLING 



A. General formalism 

We give here a brief introduction to the coupled-channel formalism for the charge ex- 
change (p, n) reaction to isobar analog states, and interested readers are referred to Ref. [jj] 
for more technical details. Let us restrict our consideration to a given isospin multiplet with 
fixed values of isospin t for the projectile and T for the target. Then, the isospin projections 
are T z = (N — Z)/2 and T z = T z — 1 for the target nucleus A and isobaric analog nucleus 
A, respectively. We further denote, in the isospin representation, state formed by adding a 
proton to A as \pA > and adding a neutron to A as \nA >. The transition matrix elements 
of the isovector part of the nucleon-nucleus optical potential ([1]) can then be obtained 
for the elastic nucleon-nucleus scattering as 

< tA\4U x {R)^-\tA >= ±^T Z U 1 (R), with r=p,n. (3) 

The + sign in the right-hand side of Eq. ([3]) pertains to incident neutron and - sign to 
incident proton. Similarly, the transition matrix element or (p, n) form factor (FF) for the 
(AT = 1) charge exchange A g . s .(p, u)Aias reaction is obtained as 

< nA\W x {R)^-\pA >= F pn (R) = jy/2T t Ui(R). (4) 

In the two-channel approximation for the charge exchange (p, n) reaction to IAS, the total 
wave function is written as 

* = \pA > Xp a{R) + \nA > x n A( R )^ ( 5 ) 

where the waves x{R) describe the relative motion of the scattering system. Then, the 
elastic (p,p) scattering and charge exchange A g , s ,(p, ti)Aias cross sections can be obtained 
from the solutions of the following coupled-channel equations [31 

[K p + U p (R) - E p ] XpA (R) = -F P n(R)x n A(R), (6) 
[K n + U n {R) - E n ] XnA (R) = -F pn (R) XpA (R). (7) 

Here K p ^ and E p ^ are the kinetic-energy operators and center-of-mass energies of the 
p + A and n + A partitions. The OP in the entrance (p + A ) and outgoing (n + A ) channels 
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are determined explicitly through the isoscalar (Uq) and isovector (Ui) parts of the nucleon 
optical potential (|T|) as 

U P (R) = U (R)-jT z Ui(R), (8) 

U n (R) = U (R) + j(T z - l)Ui(R). (9) 

In the CC calculation, both U p and U n are added by the corresponding spin-orbital potential 
as well as U p added by the Coulomb potential of the p+A system. Since the energies of isobar 
analog states are separated approximately by the Coulomb displacement energy, the (p, n) 
transition between them has a nonzero Q value. To account for this effect, the isoscalar Uq 
and isovector U\ potentials used to construct F pn [R) and U n (R) are evaluated at an effective 
incident energy of E = E^ — Q/2, midway between the energies of the incident proton and 
emergent neutron, as suggested by Satchler et al. 111]. 

Since the existing global OP parameters can be used to construct the isoscalar 

and isovector components of the proton-nucleus OP at the considered energies, it is straight- 
forward to test those parameters in the description of the A g s (p, n)AiAs reaction to isobaric 
analog states. However, more interesting structure information can be obtained when C/o(i) 



are evaluated microscopically using an appropriate folding approach 
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B. Folding model 

In our version of the single-folding model, the central nucleon-nucleus potential V 
is evaluated Hartree-Fock-type potential 

v = > + < ij\vEx\ji >], (10) 

j£A 

where vb and vex are the direct and exchange parts of the effective NN interaction between 
the incident nucleon % and nucleon j bound in the target A. The antisymmetrization of the 
nucleon-nucleus system is done by taking into account the knock-on exchange effects. To 
separate the isovector part of V which gives rise to the Lane potential, one needs to make 
explicit the isospin degrees of freedom. Namely, the following spin-isospin decomposition of 
the (energy- and density dependent) iviV interaction is used 

^d(ex) (E, p, s) = %, (EX) (E, p, s) + v® (EX) (E, p, s) (o-a- 1 ) 

+v °} EX \E, p, s)(tt') + v^ EX \E, p, s)(<t<t')(tt'), (11) 



where s is the internucleon distance and p is the nuclear density around the interacting 
nucleon pair. The contribution from the spin dependent terms (t>io and vu) in Eq. ffTTj) 
to the central nucleon-nucleus potential (1TP1) is exactly zero for a spin-saturated target like 
those considered in the present work. 

Using the explicit proton (p p ) and neutron (p n ) densities in the folding input, the nucleon- 
nucleus potential (fTO]) can be obtained [4J in terms of the isoscalar (Vis) and isovector (Viv) 
parts as 

V(E, R) = V IS (E, R) ± V W (E, R), (12) 

where the + sign pertains to incident neutron and - sign to incident proton. Each term in 
Eq. (I12p consists of the corresponding direct and exchange potentials 

Vjs(E,R) = J {[p n (r)+p p (r)]v»(E,p,s) 
+[p n {R, r) + Pp {R, r)]v**{E, p, s)j {k{E, R)s)}d\ (13) 

V lY {E,R) = J {[p n (r)-p p (r)K(E,p,s) 
+[p n (R, r) - Pp (R, r)]v**(E, p, s)j Q (k(E, R)s)}d\ (14) 

where p(r,r') is one-body density matrix of the target nucleus with p{r) = p(r, r), jo(x) is 
the zero-order spherical Bessel function, and the local momentum of relative motion k(E, R) 
is determined from 

k 2 (E, R) = ^[E c . m . - V(R) - Vc(R)}. (15) 

Here, p is the nucleon reduced mass, V(R) and Vc(R) are, respectively, the central nuclear 
and Coulomb potentials in the entrance channel (Vc = for the neutron-nucleus system). 
More details of the folding calculation of Vis and Viv can be found in Ref. [4]. 

We have further used in the folding calculation the density dependent CDM3Y6 interac- 
tion [27] which is based on the original M3Y interaction deduced from the G-matrix elements 



of the Paris NN potential 



30t ] . The densit y de pendence of the isoscalar part of the CDM3Y6 



interaction was introduced earlier in Ref. 



in numerous folding model analyses 27 



27] and its parameters have been carefully tested 
321 ] of the elastic, refractive nucleus-nucleus and 



a-nucleus scattering. Since the isovector part of the interaction can be probed in a folding 
model analysis of the charge exchange reaction only, we have developed in the present work 
an accurate procedure to parameterize the isovector density dependence of the CDM3Y6 



interaction based on the BHF results for the energy and density dependent nucleon OP 

I. The details of the new 



in nuclear matter by Jeukenne, Lejeune and Mahaux (JLM) 
treatment of the isovector density dependence of the CDM3Y6 interaction are discussed in 
Sec. HVl below. 

Given the isovector folded potential (1141) determined entirely by the neutron-proton dif- 
ference in the nuclear densities, it is necessary to have the nuclear densities determined as 
accurate as possible for a good prediction of the Lane potential. In the present work we 
have used for the considered targets the microscopic ground-state densities given by the 



Hartree-Fock-Bogoljubov approach 33j where the single particle basis includes also the con- 
tinuum states. All the results of the optical model (OM) analysis of elastic nucleon-nucleus 
scattering and CC calculation of the A gB _(p,n)AiAs reaction have been obtained with the 



CC code ECIS97 written by Raynal 



III. PREDICTION BY THE GLOBAL OPTICAL POTENTIAL 

To study the (p,n) reaction based on CC equations ©-([7]), one needs to determine the 
nucleon OP in the entrance (U p ) and outgoing (U n ) channels as accurate as possible. Since 
the elastic neutron scattering on a target being in its excited IAS cannot be measured (most 
of IAS's are either a short-lived bound state or an unbound resonance), we have determined 
U n from the isoscalar Uq and isovector U\ parts of the proton-nucleus OP evaluated at the 
effective incident energy E = _E lab — Q/2, using Eq. (J9]). The existing nucleon-nucleus global 



OP's 



0,0, Is 



8] have been carefully determined based on large experimental databases of both 
the elastic nucleon-nucleus scattering and analyzing power angular distributions, and it is 
natural to use them to construct U p for our study. The OM description of the elastic proton 
scattering from 48 Ca, 90 Zr, 120 Sn, and 208 Pb targets at incident proton energy of 40 MeV 
given by the three global proton-nucleus OP's are shown in Fig.[T]together with the measured 



data 



40 



4l| . Except some underestimation of the calculated elastic cross section in p+ 48 Ca 
case, the overall OM description of the considered elastic scattering data is reasonable. It 
should be noted that the isovector strength of the nucleon-nucleus OP is only about 2-3% 
of the total OP and its contribution to the elastic scattering cross section is too weak to 
allow us to probe the isospin dependence of the OP directly in the OM analysis of elastic 
scattering. Therefore, in a "Lane consistent" approach, the only probe of isospin dependence 
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FIG. 1: (Color online) OM description of the elastic proton scattering from 48 Ca, 90 Zr, 120 Sn, 
and 208 Pb targets at 40 MeV obtained with the global OP by Becheetti and Greenlees (BG) (J, 
by Varner et al. (CH89) [7], and by Koning and Delaroche (KD) [8]. The data were taken from 
Refs. 



of the nucleon-nucleus OP is the charge exchange v4 g . s .(p, u)Aias reaction to IAS. In such 
a quasi-elastic scattering, the charge exchange form factor (|4]) used in the CC equations 
©-0 is determined entirely by the Lane potential U\. As a result, any variation of the Ui 
strength and shape can sizably affect the calculated (p, n) cross section. Although all the 
three global OP's give about the same OM description of the elastic proton scattering as 
shown in Fig. [TJ, their descriptions of the charge exchange A g . a .(p, n)AiAs reaction are quite 
different (see Figs. [2] and [3]) . 

As discussed above, the isospin dependence of the nucleon global OP has been 
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e. m .(deg) 



FIG. 2: (Color online) CC description of the charge exchange ^4 g . s . (p, n)Au±s reaction measured 
with 48 Ca, 90 Zr, 120 Sn, and 208 Pb targets at 35 MeV obtained with Ui deduced from the global OP 
by Becheetti and Greenlees (BG) Q|, by Varner et ah (CH89) Q], and by Koning and Delaroche 



(KD) [8]. The data were taken from Ref. 



ul] 



determined from systematic OM studies of the elastic scattering of proton and neutron 
from the same target (at about the same energy), without any link to the A g s (p,n)Au^ s 
reaction. Given the absolute (p, n) cross section at least one order of magnitude weaker 
than the elastic (p, p) cross section, some discrepancy in the CC description of the (p, n) 
reaction using the Lane FF seen in Figs. [2H3] is not unexpected. From the three global 



OP's, Ui determined from the systematics by Becheetti and Greenlees (BG) js] is energy 
independent, and we found it too strong for the strength of the charge exchange form factor 
(jl]), especially, at energy of 45 MeV. Such a deficiency of the BG parameters for U\ was 
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c. m .( de 9) 



FIG. 3: (Color online) the same as Fig. [2] but for the (p, n) data measured at the proton energy of 
45 MeV. 



also found in the extensive OM analysis of elastic nucleon-nucleus scattering [7 



The 



isovector parts of both the g 
Koning and Delaroche (KD) 



obal optical potentials by Varner et al. (CH89) [7[, and by 
were found to be energy dependent and weaker than that 
given by the BG systematics. Although the KD global OP is more recent and covers a much 
wider range of energies and target masses, from the CC results shown in Figs. [2H31 one can 
see that the description of A g s (p, n)AiAs reaction by the KD global OP is slightly worse 
than that given by the CH89 global OP. A plausible reason is that the CH89 systematics was 



developed 



7J] with less constraints, based only on the elastic scattering data for A w 40 — 209 



and energies of 16 to 65 MeV (for proton) and 10 to 26 MeV (for neutron). Although 
this range of energies and target masses is narrower than that covered by the KD global 
systematics [8(, it includes the proton- nucleus systems considered in the present work. In 
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general, the Lane form factor (j4j) determined from the CH89 global OP gives a reasonable 
description of the A g _ B _(p, n)A\As cross sections measured for 120 Sn and 208 Pb targets, and 
slightly underestimates the data for 48 Ca and 90 Zr targets. As will be shown below, such a 
suppression of the calculated A g s (p, ti^ias cross sections for the two lighter targets is due 
mainly to an enhanced absorption given by the CH89 global OP. 



IV. FOLDING MODEL ANALYSIS 



A. Isovector density dependence of the CDM3Y6 interaction 



While the isoscalar density dependence of the CDM3Y6 interaction has been well tested 
in the folding model analysis 



27 



32l | of the elastic, refractive a-nucleus and nucleus-nucleus 



scattering, its isovector density dependence can be probed in the study of the charge ex- 
change reaction only. In a recent work in this direction 26| . we have used the same functional 
form for both the isoscalar and isovector density dependences and then fine tuned the scaling 
factor of the isovector part to fit the calculated (p, n) cross section to the data. Although 
we could reach good description of the (p, n) reaction under study |26j, it remains desir- 
able to have a more accurate assumption for the isovector density dependence based on 
the microscopic many-body calculation of nuclear matter. Therefore, we have developed 
in the present work a compact method to construct the isovector density dependence of 
the CDM3Y6 interaction based essentially on the BHF description of the nucleon OP in 
nuclear matter by Jeukenne, Lejeune and Mahaux 



. We recall that the isoscalar density 



dependence of the CDM3Y6 interaction was introduced 



23] as 



v^ EX) (E,p,s) =g(E)F ls (p)v^ EX \s), 
Fis(p) =C [1 + a exp(-A)P) - lop]- 



(16) 
(17) 



Parameters of the isoscalar density dependence -Fis(p) were chosen 23] to reproduce the NM 
saturation properties, with a nuclear incompressibility K pa 252 MeV, in the Hartree-Fock 
(HF) calculation of symmetric NM. These parameters as well as those corresponding to 
other K values can be found in a recent review on the nucleus-nucleus potential j^]. The 
'intrinsic' energy dependence of the isoscalar interaction is contained in the linear factor 
g(E) pa 1 — 0.0026.E, where E is the energy of incident nucleon. Given the success of the 
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TABLE I: Yukawa strengths of the central components of the M3Y-Paris interaction (|19j) 
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2.5 
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0.0 


0.0 


-7.8474 


2.6157 



parametrization ffT6]) - ffl~Tj) in numerous folding calculations, we have assumed in the present 
work a similar form for the isovector density dependence of the CDM3Y6 interaction 



(18) 



The radial shapes of the isoscalar and isovector interactions were kept unchanged, as derived 



12H from the M3Y-Paris interaction 



U 00(01) v b ) 



301 ]. in terms of three Yukawas 



Z^ r oo(oi) \y) d 
u=i UuS 



(19) 



EX 
01 



One can see from the Yukawa strengths tabulated in Table [J that the exchange terms Y c 
of the isovector interaction are much stronger than the direct terms (which is due to a 
cancellation between the even- and odd-state components). Therefore, an accurate evalua- 
tion of the exchange part of the isovector potential (jUj) is very essential in the folding model 



analysis of the (p, n) reaction. Such a domination of the exchange term in the isovector 

n 

interaction has been first noted by Love [311 ] . In our folding approach [4J the exchange parts 
of both the isoscalar (jTBl and isovector (fl4"I) proton-nucleus potentials are evaluated using 
the finite-range exchange interaction ■Wq ^ 01 ^(s) and are, therefore, more accurate than those 
given by a zero-range approximation for the exchange term. r 
Since the nucleon OP in nuclear matter can be defined 



35 



361 ] as the antisymmetrized 



matrix elements of the effective iViV interaction between the incident nucleon and those 



bound in the Fermi sea, it is given by the same h 
plane waves for the single-nucleon states 37 



artree-Fock-type potential (fTUl) . but using 
|. To determine the isovector density de- 



pendence, we have further adjusted the nucleon OP obtained with the CDM3Y6 interaction 



[m the NM limit) to reproduce the JLM density- and isospin dependent nucleon OP 



28]. 



Since the JLM potential is complex, we have used two different CDM3Y functionals to match 
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separately the real and imaginary parts of the isovector CDM3Y6 potential to those of the 
JLM potential. Namely, 



F$(E, p) = (%{Ef)[l + a^E) eM-ff(E)p) - lt{E)p], 



(20) 



so that the real (u = V) and imaginary (u = W) parts of the isovector CDM3Y6 interaction 
are determined as 



v^ EX \E,p,s) = F^E,p)v^ EX \s), 
w»r\E,P,s) = F$(E,p)v»r\s). 



(21) 
(22) 



Using Eq. (j2ip . the isovector part of the real nucleon OP in the NM limit is given explicitly 
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FIG. 4: (Color online) real part Vjy(E, p) of the isovector nucleon OP given by the isovector density 
dependence (I20p of the CDM3Y6 interaction in comparison with the JLM results 2a] at E = 35 
and 45 MeV. 



as 



V w (E,p) = F^(E,p){(p n - p p )jv + / [pj^r) - pJi{k p F r)]v**(r)j (kr)d 3 r}. (23) 
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He, * . J^r, , W - «.(-),- - AW - ,e , 8 , 0rd e r „ Be S8 e, fun , 
tion; p n and p p are the neutron and proton densities of asymmetric NM with a total density 
P = Pn + Pp and the corresponding Fermi momenta k^^ = (37T 2 p p ( n )W 3 . The momentum 
/c of the incident nucleon of mass m is determined self-consistently 37| from the nucleon 
incident energy E and real OP as 



2m 

k=^ — {E-[V ls (E,p)±V w (E,p)]}. (24) 

Here Vis(E,p) is the isoscalar part of the real nucleon OP, the (+) sign pertains to incident 
neutron and (-) sign to incident proton. Due to the self-consistent definition ( |24l) of the mo- 
mentum k, the isovector potential (1231) is obtained by an iterative procedure. After Viy(E, p) 
is determined, the isovector part Wiv(E, p) of the imaginary nucleon OP is obtained from 
the same Eq. (1231) . but with F^(E,p) replaced by F$(E,p). Our approach is to find re- 




FIG. 5: (Color online) the same as Fig. 0] but for the imaginary part Wjy(E, p) of the isovector 
nucleon OP. 

alistic parameters of the isovector density dependence fT20|) of the CDM3Y6 interaction by 
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a x 2 -fitting procedure which gives the isovector part of the nucleoli OP as close as possible 
to that of the JLM nucleon optical potential V l y' M (E, p) tabulated in Ref. [28j. To keep a 
good accuracy of this fitting procedure, instead of introducing an energy dependent scaling 
factor like g(E) in Eq. (|16|) . the density dependent parameters in Eq. (120|) have been ad- 
justed separately at each energy. As illustration, the real Viy(E, p) and imaginary Wiy(E, p) 
parts of the isovector nucleon OP at 35 and 45 MeV given by the best-fit density depen- 
dent parameters (120]) are compared with the JLM results |28j in Figs. H] and respectively. 
For each target nucleus, the parameters of complex isovector density dependence have been 
searched individually at the effective incident energy E = E p — Q/2 for the calculation of 
the (p,n) form factor (J3J) and OP in the outgoing channel In all cases, the isovector 
nucleon OP given by the best-fit parameters agree closely with the JLM results in the same 
way as shown in Figs. H] and [5] for E = 35 and 45 MeV. The numerical parameters of isovec- 
tor density dependence (I20j) at different energies E can be obtained from the authors upon 
request. For the HF calculation of nuclear matter, the isovector density dependence ( l20l 
of the CDM3Y6 interaction at energy E approaching zero has also been constructed based 
on the JLM results [39( for low energies (0 < E < 10 MeV). This set of density dependent 
parameters is used in the present work to calculate the density dependent NM symmetry 



energy S(p), defined in Eq. ([2]), by the HF method developed in Ref. [l2j explicitly for use 
with the isospin- and density dependent M3Y interaction. 

In the context of a fully microscopic OP, it is also desirable to have a realistic imaginary 
isoscalar density dependence for use in the folding calculation with the real isoscalar den- 
sity dependence ( fl6l) of the CDM3Y6 interaction. Thus, we define the imaginary isoscalar 
interaction based on the same density dependent functional ( 120]) as 

wT X \E,p,s) = F^(E,p)v^ EX \s), (25) 

then the imaginary isoscalar nucleon OP in the nuclear limit is given by 

W ls (E,p) = F™(E,p){pJ» + [ [p n j\(k F r) + P J 1 (e F r)}v^(r) Jo (kr)d 3 r}. (26) 



Here = J v® (r)d 3 r, and other involved variables are determined in the same way as 
those in Eq. (|2"3"]) . In a similar manner, the parameters of F^(E, p) have been searched 
at each energy to reproduce the JLM results tabulated in Ref. [28j. As an example, the 
isoscalar potential Wis(E, p) given by the best-fit parameters and the corresponding JLM 
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potential at E = 35 and 45 MeV are shown in Fig. [61 We must note that the imaginary OP 




FIG. 6: (Color online) imaginary part Wis(E, p) of the isoscalar nucleoli OP given by the isoscalar 
density dependent interaction ()25[) in comparison with the JLM results [28f] at E = 35 and 45 MeV. 

based on the JLM results for nuclear matter describes the absorption due to to the Pauli 
blocking effect which leads to a finite mean-free path of nucleons in nuclear medium. As 
a result, Wis(E, p) tends to have a volume shape (deep in the center and shallow at the 
surface). In general, the imaginary nucleon OP at low and medium energies has been found 
[3, s| to be best represented by a combination of volume and surface terms. The surface 
absorption is caused mainly by the collective excitations and transfer reactions which occur 
in the nuclear surface and are not related to the "volume" absorption given by W\§(E, p). 

In conclusion, we have used the HF method fl23|) - fl26|) to construct a complex isovector 
density dependence o 



in nuclear matter 



28 



the CDM3Y6 interaction based on the JLM results for the nucleon OP 
. In connection with our study, we recall that the original version of 
the JLM interaction was simply deduced from the JLM nucleon OP in nuclear matter using 
a local density approximation and smeared out in the coordinate space by a single Gaussian 



28. 



44| . For example, the real part of the original JLM interaction was constructed by this 
method from the real nucleon OP in nuclear matter V^^(E, p) as 

v 00 (oi){E,p,s) ~ exp(-— ), (27) 
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with the Gaussian range t chosen to give a good global fit to the elastic data. Since 
^is(iv) i P) a l rea dy contains strengths of both direct and exchange parts of the G-matrix, 
the nucleon-nucleus OP for finite nuclei is given by the direct folding integration (1131) only. 
Despite the simplicity, the original JLM interaction has been used quite successfully to study 



the elastic nucleon-nucleus scattering [44) as well as the (p,n) reaction to IAS [45l. 1461] . 



B. Results and discussions 

Given the new complex density dependence of CDM3Y6 interaction, the isoscalar and 
isovector parts of the nucleon-nucleus OP can be calculated explicitly by the single-folding 
approach fTl3|) - frn|) . It is natural, as the first step, to check the OM description of elastic 
nucleon-nucleus scattering at the nearby energies using the complex microscopic OP 

U(R) = N y [V IS (R) ± V W (R)} + iN w [W ls (R) ± W W (R)\, (28) 

where the (+) sign pertains to incident neutron and (-) sign to incident proton. Note that 
the imaginary part Wis(iy)(R) of the OP is given by the same folding procedure (lT5j) - ffT4"l) but 
using the imaginary parts (|22|) and (1251) of the CDM3Y6 interaction constructed separately 
at each energy. U is further added by the spin-orbital potential (and proton-nucleus OP 
added also by the Coulomb potential) taken, for simplicity, from the CH89 systematics jj]. 
The strengths iVv(w) of the complex folded OP are adjusted to the best OM fit to the elastic 
scattering data. The OM results obtained for the elastic proton scattering at 40 MeV on the 
targets under st udy are shown in Fig. [7J A good description of the measured elastic proton 



scattering data 



40 



41] can be reached after the complex folded potential is renormalized 
by Ny Rj 0.90 — 0.94 and iV w « 0.6 — 0.8. The OM results obtained for the elastic neutron 
scattering are shown in Fig. [8] where the best-fit Ny w 0.9 and Ny/ « 0.6 — 0.7. To have 
accurate distorted waves for the CC study of A g _ s _(p, n)AiAs reaction, we have also tried a 
hybrid choice for the complex OP with its real part given by the folding model and imaginary 
part by a Woods-Saxon (WS) potential based on the CH89 global systematics 

U(R) = N v [V ls (R)±V w (R)]-i[WJ(R)-4a w W s ^^}, 
where f(R) = 1/{1 + exp[(f2 - R w )/a w }. (29) 

The normalization factor Ny of the real folded potential as well as strengths of the volume 
(W v ) and surface (W s ) terms of the absorptive WS potential are fine tuned in each case to 
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FIG. 7: (Color online) OM description of the elastic proton scattering from Ca, Zr, Sn, and 
208 Pb targets at 40 MeV obtained with the complex folded OP (Hi and hybrid OP ([29]). The data 

nn 

were taken from Refs. [40I. I41I] . 



fit the elastic scattering data under study and/or reproduce in the OM calculation the total 
reaction cross section <7r measured for the considered proton-nucleus systems at 35 and 45 
MeV (43). The OM descriptions of the elastic proton and neutron scattering given by such 
a hybrid choice for the nucleon OP are shown in Figs. [7J and [8] as solid curves. One can see 
that the OM fit given by the hybrid OP is slightly improved, especially, at forward scattering 
angles. Although the difference in the OM description of elastic nucleon scattering by the 
two choices of OP is marginal as seen in Figs. [7J and El their effect on the calculated (p, n) 
cross section is surprisingly much more significant. 

After the OP for the entrance proton-nucleus channel is determined based on the OM 
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FIG. 8: (Color online) OM description of the elastic neutron scattering from 90 Zr, 120 Sn, and 208 Pb 
targets at energies of 17 to 24 MeV obtained with the complex folded OP and (|28p and hybrid OP 



29]). The data were taken from Refs. 



47|, 



48 



49fl. 



analysis of the proton elastic scattering at the nearby energies, the (complex) charge ex- 
change FF for the (p, n) transition channel is determined from the real and imaginary parts 
of the folded isovector potential ffl4"|) . evaluated at E = — Q/2, as 

2 



F pn (R) = ^/WjJ^R) = JT[N R V W (R) + iNxW w (R)). 



(30) 



Keeping the OP parameters unchanged as fixed from the OM calculation described above, 
the normalization factors N^m of the folded charge exchange FF were adjusted for the best 
fit of the calculated (p, n) cross section to the data. In this way, the folding model analysis of 
the (p, n) reaction can serve as a good probe of the isospin dependence of the effective iViV 
interaction. Since the elastic neutron scattering on a target being in its excited IAS cannot 
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be measured, the complex OP for the outgoing n + A channel has been determined from 
the complex proton OP evaluated at the effective incident energy E = Ey ah — Q/2, based 
on the isospin coupling (Q. For consistency, the complex folded OP in the n + A channel is 
renormalized by the same factors iVv(w) as those used in entrance proton-nucleus channel. 
The WS imaginary part of the hybrid OP (1291) in the outgoing n + A channel is determined 
from the CH89 global OP using the same isospin coupling (Q. The OP parameters used 
in our CC calculation of the (p, n) reaction are given in Tables [III and II I II for the complex 
folded and hybrid OP, respectively. 

We discuss now in details the CC results for the (p,n) reaction measured with 48 Ca 
target. The OM descriptions of the elastic p+ 48 Ca scattering data at 35 MeV 4QJ given by 
the complex folded OP (1281) . hybrid OP (|29|) and CH89 global OP [7| are shown in lower part 
of Fig. [9j Similar to the results at 40 MeV shown in Fig. [71 both complex folded and hybrid 
OP give a reasonable description of the measurement after their strengths were adjusted by 
the OM fit to the elastic data, with the complex folded OP slightly underestimating data 
at the forward angles. The CH89 global OP [7| sizably underestimates the elastic scattering 
data, and this is caused by a stronger absorption given b y th e CH89 imaginary OP. The 
CC descriptions of the 48 Ca g s .(p, ?t,) 48 Scias data at 35 MeV given by the unrenormalized 
folded form factor ([301 and that deduced from the isovector term of the CH89 potential using 
Eq. (j4j) are shown in upper part of Fig.[9j where the same OP's as those used in the lower part 
of Fig. [9] were taken for the entrance channel and the corresponding U n potentials evaluated 
at E = Eiah — Q/2 taken for the outgoing channel. One can see that the unrenormalized 
folded FF gives a reasonable description of the measured (p, n) cross section at large angles 
while underestimates the data points at the forward angles. From the two choices of the OP, 
the complex folded OP ([281) gives a worse fit to the (p, n) data at forward angles. Since the 
angular distribution at forward angles is more affected by the surface part of the OP and 
given the same real folded OP used in both calculations, the difference caused by the two 
OP's should be due to different surface absorptions described by the two OP's. The role of 
absorption is also seen in the CC description of the (p, n) data by the Lane FF determined 
from the CH89 parameters (denoted hereafter as CH89 form factor). Namely, the CH89 
form factor sizably underestimates the data over the whole angular range when the OP's in 
the entrance and outgoing channels are taken exactly as given by the CH89 systematics [3] • 
The CC description by the CH89 form factor improves significantly when the best-fit hybrid 
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H ^IG. 9: (Color online) upper part: CC description of the 48 Ca g . s .(p, n) 48 SciAs reaction at 35 MeV 
111 ] given by the (unrenormalized) folded (p, n) form factor (|30p and that deduced from Eq. (J3|) 
using the CH89 parameters Q]. Lower part: OM description of the elastic p+ 48 Ca scattering at 35 
MeV [40] given by the complex folded OP ([28]), hybrid OP fl29|) and CH89 global OP flj. 

OP ( |29l) is used (see Fig. [TDJ) . Therefore, the unsatisfactory description of the (p, n) data by 
the CH89 form factor shown in upper part of Fig. [9] is caused by a too absorptive imaginary 
CH89_potential (which gives ctr ~ 1124 mb compared to the measurement a^ p m 971 ± 32 



mb 43]). We have further adjusted the complex strength of the folded FF to the best x 2_ fit 
of the (p, n) data at 35 MeV 111 ], and turns out to be around 1.3, while N\ remains close 



to unity (see lower part of Fig. [TU] and Tables [TT] and IIIII) . The deficiency of the complex 
folded OP cannot be eliminated by such an adjustment of the folded FF. A consistency check 
has also been made with all the folding and CC calculations redone using the real isovector 
component of CDM3Y6 interaction increased by a factor of 1.3, and all the calculated cross 
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FIG. 10: (Color online) CC description of the 48 Ca g . s .(p, n) 48 Sc TA s reaction [lit] at 35 MeV (lower 
part) and 45 MeV (upper part) given by the renormalized folded (p, n) form factor (|30p and that 
deduced from Eq. (gj) using CH89 global OP 7], using two choices (|25]> - (|2g]l of the OP. 

sections remain about the same (as seen in the logarithmic graphing). The effect by the 
imaginary OP becomes more substantial in the CC analysis of the (p, n) data at 45 MeV 
(upper part of Fig. [10]) . While the use of the hybrid OP (|29l) results on about the same 
best-fit iV R (i) coefficients of the folded FF as those found at 35 MeV, the complex folded 
OP (128]) gives a much larger of around 1.7 and a worse description of the (p,n) data 
at large angles. Since the real OP and complex folded FF are exactly the same in both 
calculations, the difference in the calculated (p, n) cross sections (solid and dotted lines in 
Fig. [TO]) is entirely due to the difference in the imaginary OP's. For illustration, we have 
plotted radial shapes of the p+ 48 Ca optical potential at 45 MeV in Fig. [TTJ One can see 
that the real folded potential (renormalized by iVy ~ 0.9) is quite close in strength and 
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FIG. 11: (Color online) the real and imaginary parts of the complex folded OP (|28|) for p+ 48 Ca 
system at 45 MeV, scaled by factors iVyrw) given by the OM fit of the elastic data at 40 MeV, in 
comparison with the CH89 global OP taken from Ref. 7]. The adjusted imaginary CH89 potential 
has been used in the hybrid choice f)29|) of the OP. 

shape to the global CH89 real OP. The situation with the imaginary OP is quite different: 
while the imaginary folded potential has a volume shape, the imaginary CH89 potential is 
strongly surface peaked even after its strengths W v and W s were reduced by the OM fit to 
the elastic proton scattering data. The obvious reason is that the imaginary folded potential 
( I28p has been constructed based on the imaginary nucleon OP given by the BHF calculation 
of NM and is, therefore, of a "volume" nature. As a result, the imaginary folded potential 
cannot properly account for the surface absorption caused by inelastic scattering to the 
low-lying collective excitations and transfer reactions. The renormalization factor of the 
folded FF was driven to the excessive value of ~ 1.7 by the use of the imaginary folded 
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TABLE II: Renormalization coefficients iVyrw) °f the complex folded proton-nucleus OP (|28p used 
in the entrance channel. The calculated proton total reaction cross section <tr is compared with 



the data <r^ p taken from Ref. 



• -^R(i) are the renormalization coefficients of the folded FF (|30|) 



deduced from the CC fit to the (p, n) data using the OP (J2H 



Tar tret 


E 


V 


Aw 
1 v w 


u K 


u R 




Ni 




(MeV) 






(nib) 


[nib] 






48 Ca 


35 


0.933 


0.600 


969 


971 ± 32 


1.356 


0.970 




45 


0.902 


0.630 


893 


908 ± 34 


1.738 


1.054 


90 Zr 


35 


0.893 


0.731 


1341 


1316 ± 65 a 


2.133 


0.978 




45 


0.893 


0.731 


1296 


1214 ± 59 b 


2.193 


1.043 


120 Sn 


35 


0.937 


0.828 


1605 


1668 ± 59 


2.372 


0.981 




45 


0.937 


0.731 


1588 


1545 ± 38 


2.529 


0.985 


208p b 


35 


0.916 


0.747 


1877 


1974 ± 38 


2.896 


1.018 




45 


0.916 


0.747 


1963 


1979 ± 41 


2.606 


0.985 



a Total p+ 9U Zr reaction cross section measured at E = 40 MeV; at E = 49.5 MeV 



potential and not by the weakness of the isovector interaction ( 12TI) - ( 1221) . We note further 
that this subtle "absorption" effect has been established only in the CC calculation of the 
(p, n) reaction to IAS because the elastic nucleon scattering data at the nearby energies are 
still reasonably reproduced with the imaginary folded potential (see Figs. [7] and [8]). Thus, 
the distorted waves XpA and x n A g iven by the realistic hybrid OP (|2"§|) should be more 
accurate for the probe of the isovector density dependence in the CC analysis of the (p, n) 
reaction. The CC calculations using the hybrid OP (1291) give a good overall description of 
the 48 Ca g s .(p, n) 48 SciAs data at 35 and 45 MeV with the folded FF renormalized by A^ R « 1.3 
and N\ ~ 1. These calculations also give the total (p, n) cross section a pn w 10.7 and 9.0 mb 
for the 48 Ca g . s .(p, ^) 48 Scias reaction at 35 and 45 MeV, respectively, which agree well with 
the measured values [111 ]. <jp* p ~ 10.2 ± 1.1 and 8.4 ± 1.0 mb at 35 and 45 MeV, respectively. 



The results of our folding model analysis of the 90 Zr g . s .(p, n) 90 NbiAs reaction at the same 
energies are compared with the data in Fig. [121 One can see that the peak of the (p, n) cross 
section is weaker and only around half of that measured for 48 Ca g . s .(p, n) 48 SciAs reaction. 



25 



TABLE III: Parameters of the hybrid OP (|29[) used in the entrance and exit channels. Parameters 

culated 

4311 ■ ^R(l) 

are the renormalization coefficients of the folded FF (|30p deduced from the CC fit to the (p, n) 
data using the OP (1291) . 



given in boldface were kept unchanged as determined from the CH89 systematics [7|]. The ca 



proton total reaction cross section or is compared with the data cjp Xp taken from Ref. 



A 


E 

(MeV) 


Channel 


N\r 
i v V 


v v y 

(MeV) 


V V g 

(MeV) 


(fm) 


Ujyj 

(fm) 


{Tt> 
u K 


exp 

(Tr. 

U R 

( mb I 
\ 11LU / 


1 *K 


iVr 
i 


48 Ca 


35 


r> A- A 


0.925 


1.495 


5.432 


4.414 


0.69 


981 


971 ± 32 


1.265 


0.960 






n + A 


0.925 


1.495 


4.503 


4.414 


0.69 












45 


n + A 


0.900 


1.096 


5.358 


4.414 


0.69 


893 


908 ± 34 


1.279 


0.970 






n + A 


0.900 


1.096 


3.985 


4.414 


0.69 










90 Zr 


35 


n A- A 


0.913 


1.479 


6.060 


5.540 


0.69 


1330 


1316 + 65 a 


1.202 


0.969 






n + A 


0.913 


1.891 


5.267 


5.540 


0.69 












45 


r> A- A 


0.913 


2.434 


5.314 


5.540 


0.69 


1296 


1214 ± 59 b 


1.298 


1.081 






n + A 


0.913 


2.918 


4.721 


5.540 


0.69 










120 Sn 


35 


P + A 


0.937 


2.305 


7.792 


6.140 


0.69 


1637 


1668 ± 59 


1.203 


0.950 






n + A 


0.937 


1.686 


4.687 


6.140 


0.69 












45 


P + A 


0.937 


2.027 


6.529 


6.140 


0.69 


1570 


1545 ± 38 


1.225 


0.958 






n + A 


0.937 


2.653 


4.218 


6.140 


0.69 












35 


p + A 


0.901 


2.419 


8.729 


7.460 


0.69 


1964 


1974 ± 38 


1.201 


0.955 






n + A 


0.901 


1.127 


4.386 


7.460 


0.69 












45 


p + A 


0.901 


2.827 


6.334 


7.460 


0.69 


1998 


1979 ± 41 


1.150 


0.930 






n + A 


0.901 


1.871 


4.000 


7.460 


0.69 











a Total p+ yu Zr reaction cross section measured at E = 40 MeV; at E = 49.5 MeV. 

A weaker charge exchange strength also results on the total (p, n) cross section of about 
50% smaller than that obtained for 48 Ca target (see Table I in Ref. ll|). In terms of the 
isospin-flip transition (j4j), the charge exchange (p, n) strength is directly proportional to the 
neutron-proton asymmetry parameter 5 = (N — Z)/A and strength of the Lane potential 
U\. Indeed, the isovector folded potential Vpj(R) for the p+ 48 Ca system is about 30-40% 



larger than that obtained for the p+ Zr system at the surface distances and the asymmetry 
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FIG. 12: (Color online) the same as Fig. [10] but for the Zr g . s . (p, n) NbiAs reaction [11] 



parameter 5 ~ 0.17 and 0.11 for Ca and 90 Zr, respectively. A weaker charge exchange 
strength observed in the 90 Zr gs .(p, n) 90 NbiAs reaction is, therefore, well anticipated. Like 
the p+ 48 Ca system, the use of the complex folded OP fl28l) in the CC calculation with 
the folded FF gives a poor description of the (p, n) data, especially at forward angles (see 
Fig. [T2|) . even after its real strength is renormalized by Nr > 2 as determined from the x 2 
fit to the data. A strongly "volume" imaginary folded potential is also the reason for this 
disagreement. The same folded FF gives a much better fit to the (p, n) data when the hybrid 
OP ff29l) is used and its complex strengths need to be renormalized by just w 1.2 — 1.3 
and Ni w 1 which are close to those obtained for the p+ 48 Ca system (see Table [TTT]) . The 
CH89 form factor for the p+ 90 Zr system slightly underestimates the data if the OP in the 
entrance and exit channels are determined as given by the CH89 parameters. However, the 
CC description of the (p, n) data by the CH89 form factor is much better when the hybrid 
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OP (1291) is used. The CC calculation using the hybrid OP and renormalized folded FF gives 
the total (p, n) cross section a pn = 4.8 and 4.1 mb for the 90 Zr g . s .(p, n) 90 NbiAs reaction at 35 
and 45 MeV, respectively, which agree nicely with the data (cr^ p « 4.8 ± 0.5 and 4.4 ± 0.5 
mb at 35 MeV and 45 MeV, respectively) 11]. 
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FIG. 13: (Color online) the same as Fig. [TUlbut for the 120 Sn g . s .(p, re) 120 SbiAS reaction 



a. 



The CC results for the 120 Sn g s {p, n) 120 SbiAs an d 208 Pb g . s .(p, n) 208 Bi!As reactions are pre- 
sented in Figs. [13] and [TH respectively. Like the results obtained above for 48 Ca and 90 Zr 
targets, the use of the complex folded OP (128]) with a volume-shape imaginary part leads 
to a wrong shape of the calculated (p, n) cross section at forward angles (see, in particular, 
Fig. fTBl) . The CC description of the (p, n) data by both the folded FF and CH89 form factors 
is very satisfactory when the hybrid OP's (which describe well the proton elastic scattering 
at 40 MeV and measured total reaction cross section) are used for the entrance and exit 
channels. A stronger proton-nucleus Coulomb potential seems to push the main peak of the 
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(p,n) cross section to the forward angles (compare, e.g., Figs. [TOl and fl~4"|) . but the measured 
data points in the observable angular range still allow us to make an accurate conclusion on 
the complex strength of the folded (p,n) form factor (I5U|) . For the two heavy targets, the 
best CC fit by the folded FF is reached when its real and imaginary strengths are scaled 
by iVft « 1.2 and Nj « 1 which are reasonably close to those obtained for 48 Ca and 90 Zr 
targets. Although the complex folded OP (|28|) can be used to reasonably describe the elastic 
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FIG. 14: (Color online) the same as Fig. [T0l but for the 208 Pb gs .(p, n) 208 BiiAs reaction 




proton and neutron scattering on the targets under study, the volume absorption given by its 
imaginary part strongly affects the distorted waves XpA and x n A a ^ the surface an d, hence, 
leads to a poor description of the (p, n) data at forward angles and a much stronger renor- 
malization of the folded FF. In general, a fully Lane consistent and accurate description of 
both the nucleon elastic scattering and (p, n) reaction to IAS should be reached with a more 
accurate microscopic model for the imaginary OP, like that developed for the nucleon OP 
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at low energies [50|, capable to take into account explicitly coupling to the most collective 
particle-hole excitations in the target which lead to the surface absorption. Given such a 
strong impact by the absorption to the calculated (p, n) cross section, the renormalization 
factors of the folded (p, n) form factor N-^m obtained with the complex folded OP should not 
be considered as due to the weakness of the isovector density dependence of the CDM3Y6 
interaction. We must, therefore, rely on the best-fit Njim coefficients obtained with the 
hybrid OP fl29|) in deducing the strength of the isovector density dependence. Our present 
results emphasize that an accurate determination of the imaginary nucleon-nucleus OP is 
very important, especially, in the DWBA or CC studies of direct reactions measured with 
unstable nuclei when elastic scattering data are not always available. 

In connection with the present study, we note that the original version of the effective JLM 
interaction ff27|) has been used by Pakou et al. |45j and Bauge et al. [46j to study the same 



(p, n) reactions. The JLM folding-model analysis of the proton, neutron elastic scattering 
and (p, n) charge exchange reaction done in Ref. has also shown that the isovector part 
of the JLM interaction is too weak and a very strong overall renormalization of the folded 
FF by iVft = Ni ~ 2 — 2.5 is needed to account for the measured (p,n) cross sections. In 
view of our results obtained with the complex folded OP ff28l) . it is very likely that such 
large renormalization of the folded FF has been driven by the "volume" absorption of the 
JLM complex OP used in Ref. 4s|. In a more elaborate treatment of the charge exchange 
transition within the JLM model, Bauge et al. 46J have made the isospin coupling factor in 
Eq. (j4j) density dependent, i.e., ^/2T Z /A = \J[p n (r) — p p (r)]/p(r), and included it into the 
(direct) folding integral. The JLM nucleon OP obtained in such a density-dependent isospin 
coupling assumption has been thoroughly tested in the OM analysis of the proton, neutron 
elastic scattering and (p, n) reaction over a wide range of energies and target masses, and 
one can deduce from the results shown in Fig. 1 of Ref. 46j the best-fit renormalization 
coefficients of the (p, n) folded form factor « 1.5 — 1.6 and Ni w 1.3 — 1.4, in the energy 
range of 30 — 40 MeV, which are closer to our results. Despite differences in the best-fit 
renormalization coefficients of the folded FF obtained in the present work and in the JLM 



folding-model analyses 



45 



461 ] . all the results show consistently that the isovector strength 



of the JLM interaction is much too weak to account for the measured (p, n) data. Since 
the isovector term of the JLM nucleon OP has been obtained as the first-order expansion 
of the mass operator of symmetric nuclear matter perturbed by a neutron excess 5l|, a 
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weakness of the resulting JLM nucleon OP in asymmetric NM could well be expected. As 
the charge exchange reaction to IAS is quite helpful in probing the isospin dependence of 
the effective NN interaction, it would be of interest to apply similar folding model analysis 
to test the isospin dependence of the nucleon OP given by more advanced BHF calculations 
of asymmetric NM, like that by Zuo et al. 29]. 
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dependence of the nuclear symmetry energy S(p) given by the HF 



FIG. 15: (Color online) density 
calculation of asymmetric NM [12| using different isovector density dependences of the CDM3Y6 
interaction and the empirical values deduced from the CC analysis of the p( 6 He, 6 Li*)n reaction 



261 ] as well as the neutron-skin 
the text. 



2l| and HI fragmentation 



23 



24l | studies. See more discussion in 



As mentioned above, the knowledge about the isospin dependence of the in-medium NN 
interaction is of vital importance in studying the equation of state of asymmetric NM, the 
key input for the understanding of the dynamics of supernova explosion and neutron star 
formation 14], Q, [3] . We show here that the folding model analysis of the (p, n) reaction 
can be quite helpful for the determination of the nuclear symmetry energy. After the real 
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isovector density dependence (1201) of the CDM3Y6 interaction at energy approaching zero 
las been carefully parameterized on the HF level to match the microscopic BHF results 



391 ] . the density- and isospin dependent CDM3Y6 interaction is further used to calculate 



the nuclear symmetry energy S(p) using the standard HF method 12J and the results are 
shown in Fig. [Tj3 One can see that E sym = S(po) is lying somewhat lower than the empirical 
value of around 30 — 31 MeV if F^y(p) is taken as determined from the BHF results by 



JLM 



281 ] . The weakness of the JLM isovector interaction is, thus, also confirmed in our 



HF calculation of asymmetric NM. A very good agreement with the empirical -Esym values is 
reached when F^(E pa 0, p) is scaled by a factor ~ 1.3 — 1.6, which is slightly larger than 
factor for F^y(E,p) at E = 35 and 45 MeV deduced from our folding model analysis of 
the (p, n) reaction. The renormalized strength of the isovector density dependence also gives 
a good agreement with the empirical symmetry energy at the half-density point, S(p ~ 0.08 



fm~ 3 ) ^18 — 22 MeV, found from recent studies 23j, |24( of the heavy-ion fragmentation in 
the same energy region as that considered in our study of (p, n) reactions. It should be noted 
that analysis of HI fragmentation data was made based on the antisymmetrized molecular 
dynamics simulations [23[ which obtained S(p ~ 0.08 fm~ 3 ) at a finite temperature of around 
3 MeV. Therefore, this value approximately agrees with our HF result for the low-density 
part of S(p) shown in Fig. [15]only if the temperature dependence of S(p) at low NM densities 
is neglected. Finally we note that our results are also complementary to the structure studies 
which relate the E sym value to the neutron skin, a method first suggested by Brown 19|. 
If one adopts a neutron-skin AR pa 0.1 — 0.2 fm for 208 Pb then a systematics based on the 
mean- field calculations [2l| gives E syra pa 27 — 31 MeV (which is plotted as solid square in 
Fig. fl~5]) . Although the folding model analysis of the (p,n) reaction has put a constraint 
on the nuclear symmetry energy S(p) at p < po, its behavior at high densities remains 
unexplored in our study due to a simple reason that the total nuclear density of the proton- 
nucleus system never exceeds po, so that the (p, n) reaction is sensitive to the low-density 
part of the isovector interaction only. In particular, when the isovector density dependence 
Ejv(p) is taken to have the same functional form as the isoscalar density dependence Fis(p), 
and scaled by a factor of 1.1 deduced from our recent CC analysis of the p( 6 He, 6 Li*)n data 



26(| , it gives nearly the same description of the symmetry energy S(p) at p < po as the newly 



parameterized isovector density dependence (see dotted curve in Fig. [15]). The two sets of the 
isovector density dependence have, however, very different behaviors at high NM densities. 
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The S(p) curves obtained with the isovector density dependence based on the BHF results 
of asymmetric NM increase monotonically with the increasing NM density. Such a behavior 
has been recently discussed as the asy-stiff density dependence of the nuclear symmetry 
energy, while the S(p) curve given by Fjy(p) having the same functional form as -Fis(p) can 
be referred to as the asy-soft density dependence (see more discussion in Ref. [521] ). Although 
some HI collision data seem to prefer the stiff density dependence of the symmetry energy 



52j, much more studies need to be done before a definitive conclusion can be made. In 
particular, a double-folding approach to study the ( 3 He,t) or ( 13 C, 13 N) reactions exciting 
IAS might allow us to test the high density part of the isovector density dependence ( 12131 ). 
due to a higher overlap nuclear density reached during the collision and, eventually, to probe 
the nuclear symmetry energy S(p) at higher NM densities. 

V. SUMMARY 

A consistent CC analysis of the charge exchange (p, n) reactions to the isobaric analog 
states of the ground states of 48 Ca, 90 Zr, 120 Sn and 208 Pb targets at the proton incident 
energies of 35 and 45 MeV has been done using the (p, n) form factors either calculated 
microscopically in the folding model 14 or determined from the empirical WS parameters of 
the existing nucleon global OP's jf], B, 8]. 

Although the isospin dependence of the CH89 global OP [7] has been established based 
only on the OM studies of the elastic proton and neutron scattering only, it can be used 
to determine the charge exchange FF for the (p,n) transition to IAS, based the isospin 
coupling (]4]). This CH89 form factor was shown to account quite well for the (p,n) data if 
the parameters of the proton OP are fine tuned to reproduce the measured elastic proton 
scattering and total reaction cross sections ctr. 

To probe of the isospin dependence of the effective NN interaction, a complex isovector 



density dependence of the CDM3Y6 interaction 



271 ] has been carefully parameterized based 



on the density dependent JLM nucleon OP 281 and used further in the folding model analysis 



of the (p,n) reaction. Like previous studies 45|, |46| using the original JLM interaction (T2~T|) . 



the present results also show that the isovector strength of the JLM interaction is quite 
weak to account for the observed (p, n) transitions. The CC results obtained with realistic 
(semi-microscopic) nucleon OP's for the entrance and exit channels have shown that the 
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real isovector density dependence needs to be enhanced by about 20 — 30% to give a good 
description of the (p, n) reaction. 

The isovector density dependence of the CDM3Y6 interaction has also been constructed 
based on the JLM nucleon OP at energy approaching zero for further use in the HF study of 
asymmetric NM. The HF calculation using this new isovector interaction gives the nuclear 
symmetry energy S(p) close to the empirical values at p < po when the real isovector density 
dependence is scaled by a factor Nr « 1.3 — 1.6. This result confirms the weakness of the 
isovector strength of the JLM interaction found in the folding model analysis of the (p, n) 
reaction at 35 and 45 MeV. The new isovector density dependence predicts a behavior of 
S(p) at high NM densities similar to what discussed recently in the literature 52| as the 
asy-stiff density dependence of the symmetry energy. 
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